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The impurity effects on the transition temperature T c and the isotope effect are examined 
in multiband superconductors with magnetic and nonmagnetic impurities, where the effect of 
Coulomb repulsion is considered. It is shown that an internal transition of the superconducting 
state is induced by impurity doping, and that the transition is accompanied by a jump of 
the isotope-effect coefficient a. In particular, the transition is illustrated in a system with two 
electron bands. In some special cases, extended Abrikosov and Gor'kov (AG) equations for T c 
and the expressions of the isotope-effect coefficient a are obtained. Possible relevance of the 
present mechanism to the experimental results of S^RuCU is discussed. 
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The isotopc-cffcct coefficient provides important infor- 
mation on the mechanism of superconductivity. For ex- 
ample, the BCS theory based on the electron-phonon in- 
teractions is supported in nontransition-metal supercon- 
ductors, such as Hg and Zn, by their isotope-effect coef- 
ficient a, which is nearly equal to 0.5. In transition met- 
als and some compounds, large deviations from a = 0.5 
have been observed, but the deviations can be explained 
in the context of electron-phonon interactions by taking 
into account strong Coulomb repulsion, anharmonicity 
of lattice vibrations, and van Hove singularity. 1 ) 

On the other hand, in exotic superconductors, such 
as high-T c cuprates, organic superconductors, S^RuOzi, 
and heavy fermion superconductors, nonphonon mech- 
anisms of superconductivity, such as spin-fluctuation 
exchange interactions, have been examined as possi- 
ble mechanisms. In such nonphonon mechanisms, the 
isotope-effect coefficient must deviate from 0.5 markedly, 
or vanish completely. The presence of the isotope shift 
(a 0) suggests the presence of a contribution to the 
pairing interaction from electron-phonon interactions. 

Therefore, it may be difficult to derive information on 
the pairing mechanism solely from a single value of a. 
Therefore, it is useful to systematically examine corre- 
lations between a and other quantities, such as impu- 
rity concentration, transition temperature, hole concen- 
tration (for high-r c cuprates), and pressure. 

We examine the impurity effect on the isotope effect 
in this study. First, we briefly review it in single-band 
systems. In the presence of impurities, T c is reduced ex- 
cept for nonmagnetic isotropic impurity scattering and 
isotropic pairing. The ratio T c /T cQ satisfies the Abrikosov 
and Gor'kov (AG) equation 2 ) 

h <3r>=*<3 + 5Si> -*<3> *•«■>• (1) 

where £ s is a scattering rate proportional to n lmp , and 
tp(x) denotes the digamma function. We have introduced 
a function $ for convenience. 



If a slight change in the atomic mass does not affect 
the mechanisms of the pairing interactions and the pair- 
breaking effect by impurities, it is plausible to assume 
that the same equation as eq. (1) 

ln( T^^)^^2 + 2 7 rT c (M-,C s ) ) ~^2) (2) 

holds in the system with an atomic mass M' slightly 
different from M. From eqs. (1) and (2), we immediately 
obtain 

— = h-—i>'(-+—)Y\ (3) 

a L 27rTc V 2 2itT c >\ ' v ; 

which has been obtained by Carbotte, Greeson and 
Perez-Gonzalez (CGP), 3) where a = -dlnT c /dlxiM 
and cto = a(M,0). It should be noted that the relation 
between a/ao and T c /T c q, which is obtained by elimi- 
nating £s from eqs. (1) and (3), is universal in the sense 
that it does not depend on the mechanism of the pairing 
interactions, properties of impurity scattering, and even 
the pairing symmetry. 

In Sr2Ru04, however, a large deviation from the CGP 
universal relation has been observed by Mao et al. 4 ^ 
When the transition temperature is lowered by defects 
(as impurities), the coefficient a rapidly increases near 
T c w 0.94Xco from the negative value a w —0.1 near 
T c = T c0 to the positive value a » 0.2. 4 ' Below T c » 
0.94T c o, a seems to follow the CGP relation. On the other 
hand, it has been experimentally suggested that T c obeys 
eq. (I). 5 -* Since the CGP relation is derived directly from 
eq. (1) as explained above, the large deviation from it 
appears mysterious. 

In this paper, we develop a formulation of the impurity 
effect on T c and a in multiband superconductors, and 
show that an abrupt change in the isotope coefficient a 
with impurity doping is possible in some conditions, even 
when T c practically obeys an AG equation. 

The impurity effect in the multiband superconductor 
has been examined by some authors. 6 ~ 8 ' Golubov and 
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Mazin examined this subject and calculated the den- 
sity of states. 6 ^ Agterberg examined the impurity effect 
in a two-band case for Sr 2 Ru04 in the context of the 
orbital-dependent superconductivity. 8,9 ) It is known that 
Sr 2 Ru04 is a multiband superconductor. 10 ' n ) 

The appreciable isotope shifts observed in Sr 2 RuC"4 4 ) 
suggest the presence of the phonon contribution. 12 ) For 
simplicity, we ignore any possible additional nonphonon 
pairing interactions, 13 ) but the present theory is appli- 
cable whether they exist or not. 14 ) It should be noted 
that anisotropic pairing can be induced by the phonon 
interactions. 15 ~ 17 ) 

We define the impurity potential which scatters elec- 
trons from A-band to X '-band by 



uxx'(r - R») + v X x>(r - R») Si • tr(r). 



(4) 



We introduce symmetry functions 7^ (j>\\) with index a 
by which the pairing interaction Vxx'(P)P') 1S written 
in a diagonal form as 

Vxx'^v') = -Y,^ ) ^\\)9 ( xl^'\p\\'). (5) 

a 

( X) 

Here, p\\ as an argument of 7^, denotes the two- 
dimensional momentum coordinate on the Fermi sur- 
face of the A-band, where band suffixes are omitted 
from it. For example, in the isotropic case, pu can 
be replaced with the polar coordinates (0,<p), and we 

can put a = s,p x ,p y , - ■ ■ . The symmetry functions 

( x) 

7a are defined so that the orthonormal condition 
fd 2 p\\px(0,p\\h { a X \p\\h { J ) (p\\)/N x (0) = 5 aa/ is satis- 
fied, where N x (0) and px(0,P\\)d 2 p\\ denote the density 
of states of the A-band at the Fermi energy and that 
within the infinitesimal area d 2 p|| near p|| , respectively, 
f x) 

In particular, we put 7s =1. 

The gap function in the X-band is expanded by the 
symmetry functions as A x (p) = J2 a ^Xal^\p\\)- The 
impurity potentials arc also expanded as 

\uxx'(p,p')\ 2 =Y.^ X) (P\O u xxnL X 'H P {) (6) 

a 

and the same for \vxx>(p, p')\ 2 , which define «|y and 
v xx'i where the pairing interactions off-diagonal with 
respect to the suffix a arc ignored for simplicity. 

Here, we neglect the momentum dependences of j^y, 
u and v, in the direction perpendicular to the Fermi 
surface, except that we introduce the cutoff energy uj c 
(~ Wd) for the pairing interactions g^xx'- 

For an "a" -wave superconductor, the gap equation is 
obtained by the Born approximation 6 ) as 



A Xa = 2ttT^ A 



XX' 



A 



X' 



X' n'=0 



with 



A Xna = A 



Xa 



XY^Yna 



= iw « + isgn(cj„)/2rix- 



(7) 



(8) 



Here, we have put T 



T c , and defined \xx> 



gQ.Nx.Q), M XY = A XY/ 
1/tix = Y,y 1 / t ixy, and 



( Y )\ 



Axy = I/2t 2 xf, 



1 



2t~ixy 
1 



,N Y (0) u {2 $+v^S(S+l) 



,(2a) 
l XY 



,(2«: 



S(S + 1) 



(9) 



2t 2 xy 

where c a = 1 and 1/3 when "a" -wave pairing is of singlet 
and triplet pairings, respectively. We introduce matrices 
M, A, and u> whose XY elements are Mxy, A^* y , and 
w^&y, respectively, and a vector A whose A element 
is Axa- Then, we could rewrite eq. (7) as A = AAA, 
with K = ZttTYZ^o^ 1 i 1 - M)- 1 . Furthermore, if we 
define 



L = A -1 + $ 



(10) 



with <f> = LI-K, L = ln(2e 7 w c /7rT c ), and the Eulcr con- 
stant 7 = 0.57721 • • • , the gap equation can be rewritten 
in a compact form 



LA = LA. 



(11) 



The matrix $ includes all information of the impurity 
scattering. With the smallest positive eigenvalue L of the 
matrix L, T c is expressed as 



2e7 r n 
oj c exp[— L\. 

7T 



(12) 



In particular, when Axy = 1/2t 2 xy = for X ^ Y, 
the XY element of $ becomes $xy = &((' x )5xy with 

^ x 2t\x 2t2xx ^ ^ 

In the absence of the interband Cooper pair hopping 
(gxY = gxxSxy), we obtain ln(T c0 x/T c ) = <&(£x) with 
Tcox = (2e 7 /7r)(jj c exp[— 1/Xxx], which are the same as 
the AG equation except that (' x includes \/2t\xy with 
X ^ Y. 

In the general case, i.e., when Axy = 1/2t 2 xy is not 
necessarily equal to for X ^ Y, it is difficult to write 
the equation in a simple form such as that in the above. 
Therefore, we examine systems with two bands A and B 
to illustrate the mechanism of the transition. The eigen- 
values of L are obtained as 
1 



L = 



tr(L) ± \j(tr{L)f - 4det(L)j . (14) 



The smaller positive one gives T c using eq. (12). The 
explicit forms of the matrix elements of $ are obtained 
by carrying out the summation over n' as 



1 



xx 



1 



;(i+p)$(«+;(i-^H(x) 

* (,AB * C,AB 



XX 



(15) 



XX 



AB 



($(Cs) - $(C0), 



with A = B, B = A, Cab = (a - (b, and (ab = 
( A — (' B , where £ = (a and (b are the solutions to 
the equation (£ — ('a)(( ~ Cb) — AabAba = 0, that is, 

a,s = |[(a + c s )±v / (cr" 



Cb) 2 



2 LVSA ' Si;; - V ^B) rAK AB A B A\ 

the right-hand side of eq. (14) can be easily evaluated 
from eq. (10). 



Thus, 
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Now, we examine the Coulomb effects on the super- 
conductivity. In the single-band system, it is known that 
the Coulomb interaction U can be taken into account 
by replacing A with A = A — /j,*, where ^* = UN(0) 
and U = U/(l + UN(0) ]n(W/u c )). We have introduced 
an effective cutoff energy W of the Coulomb interaction, 
which is of the order of the band width. Since /i* includes 
lu c , it affects the isotope effect. 

In the multiband system, it is found that the effective 
Coulomb interaction is obtained by 

r(«) 



u XX' 



U 



XX' 



(a) T X"fj(a) 
XX"^C U X"X<> 



(16) 



X'' 



with Lq = Nx"(0)ln(W/u) c ), when the Coulomb in- 
teraction is expanded in a similar manner to that in 
cqs. (5) and (6). We have assumed that U^ x , and 
are constants within the ranges |£ p |, |£ p /| < W and Ct> c , 



respectively, while outside these ranges, U 



(a) 
XX' 



and 

U xx , = 0, where £ p denotes the single-particle energy 
measured from the Fermi energy. 
In two-band systems, we obtain 

1 



U 



U 



(a) 
XX 

(a) 
XX 



D 
1 

D 



(1 + L*U£ X )U£ X 



L c U XX U 



(a) 
XX 



(1 



L cUxx)Uj["x 



T X TT (a) TT {a) 
c u XX U XX 



(17) 



rA T r(a) 



)(1 



with D = (1 

The coupling constants A^i-, should be replaced with 



^XX' 



XX' 



^XX' 



^xx 
where 



#2) 



tAtA^tBtt 



(a) 
BB • 



^XX' 



= U^l,N X '(0). The 



matrix L is redefined by L = (A — p,*) + <& with a 
matrix //* whose XY element is H* XY - 

The transition temperatures and the isotope-effect co- 
efficients can be explicitly calculated from the set of 
equations obtained above. However, for simplicity and 
in order to clarify a physical picture, we concentrate 
ourselves on a symmetric case in which two bands are 
equal, namely, Na(0) — Nb(0) = N(0), Jaa = Jbb and 
Jab = /ba, where / = g( a \ u^ a \ and f/K In this 
case, the gap functions A^ with Aau = ±Ab q are the 
solutions of eq. (11). The eigenvalues are 

1 



L(±) = A 



+ <t>(C (±) ), 



(18) 



(±) 



where A (±) = X AA ± \ A b = \±) - M( ±) with A (±) = 
>^aa±^ab and ^* ±) = /j,* AA ±[i AB , and ( (±) = C' T ^ab 
with C = C' A = C'b- From cq. (17), we have 



{U A a i±U A a l)N{Q) 



l + (U A a }±U^)N(0)ln(W/^ 



Therefore, we finally obtain 



In- 



with T£> = (2e 7 /7r)w c exp[-l/A (±) ] and 



(19) 



(20) 



->(±) 

C 

C ( ± ) = (: 



1 



1 



-) + (; 



1 



1 



-)• (21) 



v 2tiaa 2t 2 aa ^TiAB 2t 2 ab 
From eq. (20), the isotope-effect coefficients for T c ^ 



E-h" 




Fig. 1. Impurity concentration dependence of the transition tem- 
perature T c . The parameter £' is proportional to the impurity 
concentration «i mp . The solid line shows the resulting transition 
temperature for each f . The dotted line shows lower transition 
temperatures which are actually suppressed. 




Fig. 2. Isotope-effect coefficient as a function of transition tem- 
perature T c . The solid circles are the experimental data by Mao 
et aZ. 4 > 



are obtained as 



a, 



with 



(±) 
o 

„(±) 



C (±) 



2:r7: (±)^2 + 2^' " 



(22) 



|[1 - (M(±)A(±)) 2 ]. In general, 



„(±) 



largely different due to the difference in M(±)- The physi- 
cal origin of these differences is that the Coulomb energy 
U A a A is enhanced or reduced by U A °g depending on the 
signs ±, as explicitly shown in eq. (19). 

For the s-wave pairing and the s-wave impurity scat- 
tering, C (+) = 0, while C ( ~ } = 1/tiab- Thus, T c of the 
state with A As = Ag s is not suppressed by impurity 
doping at all, while that of the state with Aa s = — Ab s 
is strongly suppressed, which is consistent with the re- 
sults by Golubov and Mazin 6 ) and Arseev et al. 7 ^ When 
T c ( + ) < T c (~\ the alternation of the solution and thus 
the jump of a occurs. 

For anisotropic pairing and isotropic impurity scatter- 
ing, since 1/2t2xy — 0, we have = ■ Hence, the 
alternation of the solution does not occur. 

Even for anisotropic pairing, it is possible that Aab = 
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l/2r 2y iB 7^ 0, when u A g ^ or v A g ^ 0. It is found in 
cq. (21) that when Xab^-ab < 0, one of T c ^ with lower 
T c o decreases more gradually with impurity doping than 
the other. In this case, the gap function may alternate 
from Aau = ±Ag Q to Aa u = tAbo at some impu- 
rity concentration. In the transition from A Aa = ±Ab q 
to A Ac — tAbc, the isotope-effect coefficient a jumps 
because of the difference in mentioned above. 

In Figs. 1 and 2, we show an example. To reproduce the 
experimental data of Sr 2 RuC>4, 4 ) we put = —0.1, 

a { - ] = 0.2, T c < +) = 1.5K, T c < _) = 0.98T c < +) , and 
Aab/C = -0-2. 18) The values of T c { ) ±) and are re- 
produced, for example, by uj c = 410K, 19 ) Xaa = 0.33678, 
Xab = 0.028509, fj,* AA = 0.16276, and fj,* AB = 0.028203. 
In this case, since Xab > 0, we have > tJ \ In 

Fig. 1, the curves of T c (±) cross at (' /T c0 « 0.063 and 
T c /T c0 w 0.94. Since the fold of the curve of T c is slight 
for the present parameters, it would be hardly detected 
experimentally. At the crossing point, the superconduct- 
ing state alternates from the state with A Aa = Ab q to 
that with Aaoc — -Ag a . In Fig. 2, it is shown that in this 
internal transition, a jumps. The present theory could 
reasonably reproduce the experimental result within the 
error bar of the data. In experiments, possible inhomo- 
geneity of the samples might smear the transition. 

The internal transition and the jump of a could occur 
in more general situations than in the above example. 
The eigenstates of the linearized gap equation, eq. (11), 
are subject to different impurity and Coulomb effects as 
shown in eqs. (19) and (22) in the example. Since the 
eigenstate with the highest T c occurs, the superconduct- 
ing state alternates from one eigenstate to another with 
impurity doping for appropriate parameters. In the tran- 
sition, the isotope-effect coefficient a jumps because the 
effective Coulomb parameter \i* changes. 

Lastly, we discuss the application to Sr 2 RuC>4. The 
Sr 2 Ru04 compound has three electron bands, called a, 
/3, and 7, which have separate Fermi surfaces. 10, n ) There 
are some experimental results to support triplet pair- 
ing, 20 ) but the momentum dependence of the gap func- 
tion is controversial. 2 1_23 ) As argued above, the present 
mechanism holds, as long as the parameters satisfy the 
condition for the alternation of the eigenstates, indepen- 
dently of the number of electron bands and the momen- 
tum dependence of the gap function, whether additional 
nonphonon pairing interactions exist or not. 14 ) However, 
for the alternation at T c w 0.94T c0 to be quantitatively 
reproduced, it is necessary that two of the eigenstates 
have very close T c0 . It could occur as a result of a com- 
bination of three electron bands. The split of the critical 
field curve observed for parallel fields 24 ' may suggest an 
existence of the hiden eigenstate slightly below T c0 at the 
zero field. 

We may consider another possibility for Sr 2 Ru04. The 
internal degrees of freedom play an essential role in the 
present mechanism. In the present model, they originate 
from the multiband nature, but they may originate from 
the anisotropic gap structure of triplet pairing supercon- 



ductivity. If we apply the above calculation to a model 
with a pair of the eigenstates with slightly different T c0 's 
by replacing the band suffixes A and B with the suffix a 
to express the momentum dependence, the same figures 
as Figs. 1 and 2 are reproduced for appropriate param- 
eter values. The explicit calculation will be presented in 
a separate paper. 

In conclusion, we have examined the impurity effect 
in multiband superconductors, and obtained the expres- 
sions of T c and a. It has been found that an internal 
transition of the superconducting state is induced by im- 
purity doping in some condition, and the transition is ac- 
companied by a jump of the isotope-effect coefficient a. 
The T c dependence of a deviates from the CGP universal 
relation due to the jump. It is possible under appropri- 
ate parameters that T c appears to obey the standard AG 
equation simultaneously with a large deviation from the 
CGP relation. 

The author wishes to thank Y. Maeno and Z.Q. Mao 
for useful discussions and for sending their experimental 
data to him. 
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